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The  Inverse  Scattering  Problem 
for  Acoustic  Waves* 

by 

David  Colton 

All  important  decisions  must  be  made  on  the  basis 
of  insufficient  data  -  If  You  Meet  the  Buddha  on 
the  Road,  Kill  Himl  by  Sheldon  Kopp 

I.  Introduction. 

In  this  paper  we  shall  survey  recent  progress  and  discuss 
open  problems  connected  with  the  inverse  scattering  problem  for 
acoustic  waves.  However  before  proceeding  it  is  first  necessary 
to  be  more  precise  on  what  we  mean  by  "the  inverse  scattering 
problem"  since  this  phrase  has  been  used  to  describe  a  large 
variety  of  problems  concerned  with  target  identification  in  acous¬ 
tic  wave  propagation.  ^We  first  make  a  distinction  between 
"scattering*  and  "diffraction*  and  note  that  the  latter  is  basically 
a  high  frequency  phenomena  whereas  the  former  is  more  accurately 
applied  to  low  and  intermediate  values  of  the  frequency.  Hence-  in 
this  paper  we  shall  not  discuss  any  of  the  important  new  results 
on  the' "inverse  diffraction  problem^  but  instead  refer  the  reader 
to  the  recent  paper  of  Brian  Sleeman  for  a  survey  of  these  results 
({151).  We  shall  further  restrict  our  attention  to  the  scattering 
of  a  plane  time  harmonic  wave  by  a  fixed  bounded  obstacle  situated 
in  a  homogeneous  medium  and  in  particular  to  determine  information 

*  This  research  was  supported  in  part  by  AFOSR  Grant  79-0085. 
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^about  this  obstacle  from  a  knowledge  of  the  asymptotic  behavior 
of  the  scattered  wave.  Hence,  we  are  excluding  such  topics  as 
the  scattering  of  waves  by  moving  obstacles,  the  determination 
of  the  speed  of  sound  in  a  non-homogeneous  medium,  and  the  loca¬ 
tion  of  equivalent  sources.  Finally,  we  shall  only  be  concerned 
with  determining  two  basic  properties  of  the  scattering  obstacle, 
viz.  its  shape  and/or  its  surface  impedance.  We  note  that 
although  the  problem  of  determining  the  shape  of  an  obstacle  from 
far  field  measurements  has  long  been  recognized  as  a  basic  problem 
in  a  variety  of  areas  of  technology  such  as  radar,  sonar,  and 
tomography,  the  inverse  impedance  problem  has  received  less  atten¬ 
tion.  Nevertheless  the  inverse  impedance  problem  is  of  basic 
importance  in  many  applications  since  it  gives  information  on  the 
material  composition  of  the  unknown  scattering  object,  e.g.  in 
the  case  of  sonar  it  can  help  answer  the  auestion  of  whether  or 
not  the  scattering  obstacle  is  a  whale  or  a  submarine. 

The  inverse  scattering  problem,  as  defined  above,  is 
particularly  difficult  to  solve  for  two  reasons:  it  is 
(1)  nonlinear  and  (2)  improperly  posed.  Of  these  two  reasons 
it  is  perhaps  the  latter  that  creates  the  most  difficulty.  In¬ 
deed  we  shall  see  shortly  that  for  a  given  measured  far  field 
pattern  in  general  no  solution  exists  to  the  inverse  scattering 
problem,  and  if  a  solution  does  exist  it  does  not  depend  con¬ 
tinuously  on  the  measured  data.  Hence  before  we  can  begin  to 
construct  a  solution  to  the  inverse  scattering  problem  we  must 
answer  the  question  of  what  we  mean  by  a  "solution".  At  this 
point  it  is  worthwhile  recalling  the  remark  of  Lanczos:  "A  lack 
of  information  cannot  be  remedied  by  any  mathematical  trickery". 
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Hence  in  order  to  determine  what  we  mean  by  a  solution  it  is 

necessary  to  introduce  "nonstandard”  information  that  reflects 

the  physical  situation  we  are  trying  to  model.  Having  resolved 

the  question  of  what  is  meant  by  a  solution,  we  then  have  to 

actually  construct  this  solution,  and  this  is  complicated  not 

only  by  the  fact  that  the  problem  is  nonlinear,  but  also  the 

fact  that  the  above  mentioned  "nonstandard"  information  has  been 

incorporated  into  the  mathematical  model. 

We  shall  now  give  a  mathematical  formulation  of  the  inverse 

scattering  problem  and  outline  the  specific  topics  we  want  to 

2 

discuss  in  this  paper.  We  shall  formulate  our  problem  in  ft  and 

state,  when  appropriate,  the  necessary  modifications  that  are 

needed  to  consider  the  full  three  dimensional  problem.  Let  D  be 

2 

a  bounded,  simply  connected  domain  in  ft  with  smooth  boundary  3D 
and  unit  outward  normal  v.  If  we  let  X  *  X(x)  ^>0  denote  the 
continuous  surface  impedance  of  3D  and  k>0  the  wave  number,  then 
the  impedance  boundary  value  problem  for  acoustic  waves  can  be 
mathematically  formulated  as  the  problem  of  finding  a  function 
ueC2(ft2\D) nC1(ft2\D)  such  that 


u  *  u1  + 

us 

in 

*2\d 

(1.1a) 

*2 

u  +  k2u  ■ 

0 

in 

ft2\& 

(1.1b) 

3u 

Jv 

+  ikXu  ■ 

0 

on 

3D 

(1.10 

lim 

rl/2  (j£ 

-iku8) 

-  0 

(1. Id) 

r-**» 
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where  the  "incoming  wave"  u1  is  an  entire  solution  of  the 
Helmholtz  equation  (1.1b)  and  the  "scattered  wave"  us  satisfies 
the  radiation  condition  (l.ld)  uniformly  with  respect  to  0  where 
(r,9)  are  polar  coordinates.  We  shall  also  consider  the  case 
when  A  is  infinite,  i.e.  the  boundary  condition  (1.1c)  becomes 

u  =  0  on  3D.  (1.1c') 

In  this  case  we  shall  refer  to  the  scattering  obstacle  as  being 
"soft".  The  existence  and  uniqueness  of  a  solution  to  (1.1a)- 
(1. Id)  and  (1.1a),  (1.1b),  (1.1c'),  (l.ld)  is  well  known 
(c.f.  [10]).  We  shall  see  shortly  that  if  u  is  a  solution  of 
(l.la)-(l.ld)  (or  (1.1a),  (1.1b),  (1.1c'),  (l.ld)  then  us  has 
the  asymptotic  behavior 

us(r,0)  «  elkr  r“1/2F(0;k)+O(r“3/2)  (1.2) 

where  F  is  known  as  the  far  field  pattern  corresponding  to  the 
given  incoming  wave  u1.  The  inverse  scattering  problems  we  shall 
discuss  in  this  paper  are  (1)  given  a  knowledge  of  u1,  F  and  D, 
find  X,  or  (2)  given  a  knowledge  of  u1,  F  and  the  boundary  condi- 

V 

tion  (1.1c  ) ,  find  D.  We  shall  be  more  precise  on  what  we  mean 
by  a  "knowledge"  of  F  shortly,  but  as  mentioned  above  we  in  general 
only  know  F  from  measurements  which  are  by  definition  inexact  and 
this  fact  makes  both  of  the  above  inverse  scattering  problems 
improperly  posed.  The  basic  fact  we  begin  with  is  that  the 
existence  of  a  unique  solution  to  the  direct  scattering  problem 
(1. la) - (l.ld)  or  (1.1a),  (1.1b),  (1.1c'),  (l.ld)  by  the  method 
of  integral  equations  (c.f.  [10])  defines  a  (nonlinear)  mapping 
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T  from  D  or  A  to  P.  Hence  from  an  abstract  point  of  view  we  cam 
state  our  tasks  as  follows: 

(1)  Determine  the  range  of  T  (denoted  by  R(T) )  as  a  subset  of 
L2  [0 ,2ir]  . 

(2)  Establish  the  existence  of  on  R (T) ,  i.e.  show  the 
uniqueness  of  the  solution  to  the  inverse  scattering 
problem.. 

~-l 

(3)  Determine  a  subset  XcR(t)  and  an  operator  T  defined  on 

2  i 

L  [0,2tt]  such  that  T  =  T  on  X  and  T  is  continuous 
2 

on  L  [0,2ff],  i.e.  stabilize  the  inverse  scattering  problem 
(In  order  to  do  this  it  is  necessary  to  assume  a  priori 
information  about  D  or  A) . 

—  1 

(4)  Give  a  constructive  method  for  determining  Tx,xeX. 

We  shall  examine  what  is  known  about  the  above  four 
problems  in  the  following  sections  and  in  addition  give  directions 
that  should  be  taken  if  further  progress  is  to  be  achieved. 

II.  The  Mapping  T  and  its  Range. 

As  mentioned  in  the  Introduction  we  can  obtain  the  solution 
of  (1.1a) -(1. Id)  (or  the  corresponding  problem  with  Dirichlet 

I 

boundary  condition  (1.1c  ))  by  the  method  of  integral  equations. 

In  particular  we  can  represent  the  solution  us  of  (1. la) - (1. Id) 
in  the  form  of  a  modified  single  layer  potential 

us(x)  «  i  |  $ (y)G(x,y) ds (y)  ;  X€*2\U  (2.1) 

3D 
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where  G  is  an  appropriately  chosen  fundamental  solution  ([22]) 
and  $  is  the  unique  solution  of  an  integral  equation  of  the  form 

i 


+  Kj#  +  AK24>  =  +  Au1 


(2.2) 


where  and  K2  are  compact  integral  operators  that  are  indepen¬ 
dent  of  A  and  defined  on  C(3D).  The  solution  of  the  corresponding 
Dirichlet  problem  can  be  represented  in  the  form  of  a  double 
layer  potential  ([22]) 

us(x)  =  1  |  tp(y)  G  (x,y)  ds  (y)  ;  xefc2\I5  (2.3) 

3D 

where  is  the  unique  soluton  of  an  integral  equation  of  the  form 

+  K3<j>  =  u1  (2.4) 

where  is  a  compact  integral  operator  defined  on  COD)  .  If  in 
(2.1),  (2.2)  or  (2.3),  (2.4)  we  now  let  r=  |x|  tend  to  infinity 
and  use  the  asymptotic  behavior  of  G(x,y)  we  obtain  the  relation¬ 
ship  (1.2)  and  the  mapping  T:A-*-P  or  T:3D-*-F.  In  particular  such 
a  calculation  establishes  the  validity  of  the  following  theorem: 


Theorem  1:  The  far  field  pattern  is  an  entire  function  of  8  and 
a  continuous  function  of  k  for  k>0. 


The  above  analysis  shows  both  the  nonlinear  nature  of  the  operator 

2 

T  as  well  as  the  fact  that  R(T)/L  [0,2ir].  Furthermore  since  it  is 
not  possible  to  determine  the  analyticity  of  a  function  from 
inexact  measurements.  Theorem  1  implies  that  in  general  for  a  given 
measured  far  field  pattern  no  solution  exists  to  the  inverse 
scattering  problem  unless  further  assumptions  are  made. 


\ 
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From  (2.2)  or  (2.4)  it  is  clear  that  the  ODerator  T  depends 
on  u1,  i.e.  T=T(ux).  We  can  therefore  pose  the  following  question: 
If  u1  ranges  over  all  entire  solutions  of  the  Helmholtz  equations 
where  A  or  D  is  kept  fixed,  is  the  corresponding  set  of  far  field 
patterns  dense  in  ir[0,2-n-],  or  more  concisely,  does  R(f)=L  [  0 , 2-rr  ]  . 
The  following  example  shows  that  this  is  not  true  in  general. 

Example :  Consider  problem  (1.1a),  (1.1b),  (1.1c  ),  (l.ld)  when  D 
is  the  unit  disk.  Then  since  u1  is  an  entire  solution  of  the 
Helmholtz  equation  we  can  expand  u1  in  the  form 

u*(r,8)  =  I  J  (kr)  [a„cos  n6  +  b„  sin  n9]  (2.5) 

n=0  n  n  n 


where  J  denotes  Bessel's  function  and  the  series  (2.5)  is  uniformly 
n 

2 

convergent  on  any  compact  subset  of  R  .  Then  for  r>l  we  can  expand 
us  in  the  uniformly  convergent  series 


S  00  (1) 

us  (r,0)  =  -  Z  h'-u  (kr)  — - 

n=0  n  „<1) 


[a  cos  n0  +  b„  sin  n0]  (2.6) 

n  n 


H^A/ (k) 


where  denotes  Hankel's  function  of  the  first  kind  and  from 
n 

the  asymptotic  behavior  of  Hankel's  function  we  have  that  the  far 


field  pattern  for  us  is  given  by 
F(0  ;k)  *  e*"/  ‘  Z 


_•  _  ,,  m -  00  (-i)  J_  (k) 

'  *=-r-  Z  —r -r-\ - [a„  cos  n0  +  b  sin  n0](2.7) 

TTk  _ _ n  tr(l)  /.-v  n  n 


n*0  h^a' (k) 


2 

If  kQ  is  an  eigenvalue  of  the  interior  Dirichlet  problem  then 

for  some  integer  n=n„  and  hence  in  this  case  F(8;k  ) 
no  o  o 

is  orthogonal  to  cos  nQ0  and  sin  nQ0  for  all  incident  fields 

u*.  Hence  the  class  of  far  field  patterns  for  such  values  of  k 

2 

is  not  dense  in  L  [0,2ir], 
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It  is  an  open  question  to  determine  if  similar  examples 
are  valid  for  arbitrary  domains  D.  If  such  examples  exist  this 
would  establish  an  interesting  relationship  between  the  far  field 
patterns  of  exterior  boundary  value  problems  for  the  Helmholtz 
equation  and  the  (interior)  eigenvalue  problem  for  Laplace's 
equation. 

The  validity  of  the  above  example  is  based  on  the  fact 
that  the  set 


J  (kr)  cos  n0 
n  ;  n 

J  (kr)  sin  n6 
n 


0,1,2, 


(2.8 


2 

is  incomplete  in  L  (3D)  if  k  is  eigenvalue  of  the  interior 
Dirichlet  problem.  However  it  can  be  shown  ( [5]  ;  a  simpler 
proof  can  be  based  on  the  ideas  of  [13])  that  the  set 


(■jy  +  ikA)  J„  (kr)  cos  n9 


n 


(yy  +  ikA)  JR(kr)  sin  n9 


n=0,l,2, . . . 


(2.9) 


2 

is  complete  in  L  (3D)  for  arbitrary  bounded  domains  provided 
0<A<°°  where  A  is  a  constant.  This  fact  implies  the  following 
theorem  ([5]): 


Theorem  2 ;  Let  A  be  a  constant  such  that  0<A<°°.  Then  R(T)  = 
L2[0,2tt]  . 

A  wealth  of  questions  concerning  R(T)  remain  unanswered 
and  this  provides  a  major  mathematical  challenge  of  basic  impor¬ 
tance  to  the  inverse  scattering  problem.  We  have  already  men¬ 
tioned  one  of  these  questions  in  connection  with  the  Dirichlet 
problem.  Another  is  the  following:  Determine  the  compact  set  of 
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2 

far  field  patterns  in  L  [0,2 tt)  that  corresponds  to  a  fixed  incident 
field  and  wave  number,  but  with  (constant)  surface  impedance  lying 
in  a  given  compact  subset  of  the  positive  real  axis.  As  opposed 
to  the  problem  answered  by  Theorem  2,  this  problem  is  complicated 
by  the  fact  that  it  is  nonlinear.  A  similar  problem  can  be  posed 
for  A  fixed  but  3D  lying  in  a  compact  set  of  closed  surfaces. 

III.  The  Existence  of  T-1. 

As  we  have  already  mentioned  in  the  Introduction,  the 
existence  of  T-"*-  is  equivalent  to  establishing  the  uniqueness  of 
solutions  to  the  inverse  scattering  problem.  We  first  consider 
the  inverse  scattering  problem  of  determining  D  from  a  knowledge 
of  u1  and  F  where  we  assume  that  the  boundary  data  is  given  by 
(1.1c’)  ([17]). 

Theorem  3:  D  is  uniquely  determined  by  a  knowledge  of  the  far 
field  pattern  F  for  all  angles  0  and  k  on  any  interval  of  the 
positive  real  axis. 

Proof :  Suppose  there  existed  two  obstacles  D.^  and  D2  having  the 

same  far  field  pattern  F.  Consider  first  the  case  when  and 
I?2  are  disjoint.  Then  from  the  analyticity  of  solutions  to  the 
Helmholtz  equation,  F  uniquely  determines  us  outside  a  disk 
containing  and  D2  in  its  interior  and  we  can  conclude  by  analy¬ 
tic  continuation  that  us  is  an  entire  solution  of  the  Helmholtz 

equation  satisfying  the  radiation  condition.  But  this  implies 

s  *  i 

that  u  =0  ([15]),  i.e.  from  (1.1c  )  we  have  u  =0  on  3D^.  Hence 

u1  is  an  eigenfunction  of  the  Laplacian  for  an  interval  of  k 

values  and  this  is  a  contradiction  since  the  set  of  eigenvalues 

for  the  Laplacian  is  discrete.  A  similar  argument  holds  when 

vf 
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D=D, nDj  is  nonempty,  in  particular  by  analytic  continuation  one 
can  show  that  the  Laplacian  defined  in  D^\D  or  D2 \D  has  a  contin¬ 
uous  spectrum,  which  is  a  contradiction. 

We  now  consider  the  inverse  scattering  problem  of  deter¬ 
mining  X  from  a  knowledge  of  u1 ,  F  and  D,  i.e.  the  boundary  value 
problem  is  given  by  (1. la) - (1. Id) .  The  following  theorem  is  due 
to  Colton  and  Kirsch  ([8]). 

Theorem  4 :  X  is  uniquely  determined  by  a  knowledge  of  the  far 
field  pattern  F  for  all  angles  0  and  any  fixed  positive  value  of 
the  wave  number  k. 

Proof:  Let  X-^  and  X2  be  two  solutions  of  the  inverse  scattering 

problem  corresponding  to  the  same  far  field  pattern  F.  Then  by 
analytic  continuation  we  can  conclude  that  Uj=u2  and  9ul/3v  = 
9u2/3v  on  3D  where  u^  and  u2  are  the  solutions  of  ( 1. la) - ( 1. Id) 
corresponding  to  X^  and  X2  respectively.  Then  from  the  boundary 
condition  (1.1c)  we  have 

(Xi-X2)Ui  *  0  on  3D.  (3.1) 

We  now  note  that  u^  is  not  identically  zero  in  any  neighborhood 

Sc3D  since  if  this  were  the  case  (1.1c)  implies  that  9ul/3v=0  in 

2 

S  and  hence  by  Holmgren's  uniqueness  theorem  ([3])  u^rO  in  2  \D. 
But  this  is  a  contradiction  since  u^u^+u^  and  u®  satisfies  the 
radiation  condition  but  u^  does  not.  Therefore  if  xe3D  there 
exists  a  seouence  of  points  xn-*-x  such  that  u^(xn)^0.  From  (3.1) 
we  have  that  X^  (*n)  **2  ^xn^  and  since  and  *2  are  assumed  to 
continuous  we  have  X^(x)«X2(x).  Since  x  is  an  arbitrary  point  on 
3D  this  completes  the  proof  of  the  theorem. 


i 
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We  note  that  from  Theorem  1  and  the  identity  theorem  for 
analytic  functions,  it  suffices  in  Theorem  3  and  Theorem  4  to 
assume  that  F  is  only  known  for  an  interval  of  6  values  contained 
in  [ 0 , 2tt ]  instead  of  for  all  values  of  0  in  [0,2ir]. 

An  open  problem  of  considerable  interest  is  to  determine 
the  validity  of  the  above  theorems  if  instead  of  the  far  field 
pattern  we  are  given  the  scattering  cross  section  a  defined  by 

/2tt  9 

a  (k  ;a)  =  |F(0jk,a)  Td0  (3.2) 

J  o 

where  F(0;k,a)  is  the  far  field  pattern  corresponding  to  the 
incoming  plane  wave  u1=exp (ikx . a) . 

'■-I 

IV.  The  Operator  %  « 

As  pointed  out  in  the  Introduction,  the  inverse  scattering 
problems  we  are  considering  in  this  paper  ? re  improperly  posed  in 
the  sense  that  in  general  no  solution  exists  and  if  a  solution 
does  exist  it  does  not  depend  continuously  on  the  initial  data. 
Hence  a  major  task  in  a  satisfactory  treatment  of  these  problems 
is  to  restore  stability.  We  shall  accomplish  this  by  assuming 
extra  a  priori  information  is  available  concerning  the  unknown 
impedance  or  scattering  obstacle  such  that  it  is  possible  to  con¬ 
clude  that  the  exact  far  field  pattern  lies  in  a  compact  set 

XcR(T) .  We  shall  then  show  that  it  is  possible  to  define  an 

''-l  2  A-1  -1  "'-l 

operator  T  defined  on  L  [0 , 2tt 3  such  that  T  =T  on  X  and  T 

2 

is  continuous  on  L  [0,2ir],  thus  restoring  the  continuous  dependence 


of  \  or  D  on  the  far  field  pattern.  At  this  point  we  would  like 
to  mention  that  the  term  "a  priori  information"  is  not  meant  to 
imply  that  such  information  is  impossible  to  obtain  from  experimental 
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data,  but  only  that  this  (in  general  nonstandard)  information  is 
sufficient  to  stabilize  the  problem  under  consideration. 

We  first  consider  the  problem  of  determining  the  impedance 
A  from  a  knowledge  of  the  far  field  pattern  F  where  the  incident 
wave  u1  and  the  scattering  obstacle  are  assumed  known  ([8]).  By 
using  the  compactness  of  the  operators  and  Kj  in  (22)  it  is 
possible  to  show  that  the  mapping  T:A-*-F  is  continuous  and  in  parti¬ 
cular  if  <j>  is  the  solution  of  the  integral  equation  (2.2)  we  can 

write  F=F ( 9 ; <j> , A )  where  F  is  a  continuous  function  of  d,p  and  A. 

2 

Now  let  FeL  [0,2ir]  be  the  measured  far  field  pattern.  Then  we 
can  reformulate  our  inverse  scattering  problem  as  an  optimization 
problem,  denoted  by  P^,  as  follows:  Minimize 

/■2ir  - 

Cf  ( <J> ,  A)  =  1  F  (9  ?<J> ,  A)  -f  (0)  |  d©  (4.1) 

o 

subject  to  <j>  being  a  solution  of  (2.2)  and  AeU  where  U  is  a 
"control  set"  to  be  specified  shortly.  We  first  note  the  following 
facts : 

(1)  The  set  u  contains  the  a  priori  information  we  are 
assuming  about  A. 

* 

(2)  If  A  eU  is  a  solution  of  the  inverse  scattering  problem 

* 

for  a  given  far  field  pattern  f  and  p  is  the  corresponding 

★  *  *  ★ 
density,  then  ($  , A  )  is  a  solution  of  Pf  since  Cf(<j>  ,A  J^O. 

(3)  If  (p* , A*)  is  a  solution  of  Pf  and  Cf ($*, A*)=0  then  A 

is  a  solution  of  the  inverse  scattering  problem.  If 
*  * 

Cf(<P  r  A  )  >0  then  the  inverse  scattering  problem  is  not 
*  * 

solvable  for  A  eU ,  but  A  is  a  best  approximation  in  the 
sense  that  |  |F-f  j  I L2 [0 r 2tt] 


is  minimal. 
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We  now  define  our  control  set  U  in  such  a  way  that 
has  a  solution.  In  particular  let 

U  =  UeC+OD)  :  |A(x)  I  <M1,  1  A(x)-X(y)  |<_M2|x-y!a} 

where  C+(3D)  is  the  cone  in  COD)  consisting  of  all  continuous 
functions  A  defined  on  3D  such  that  A >0 ,  and  ,  M2  and  a,0<o<^l, 
are  fixed  constants.  The  Arzela-Ascoli  Theorem  implies  that  U 
is  compact  in  COD)  and  since  the  functional  Cf  is  continuous 
we  have  the  following  theorem: 

2 

Theorem  5 :  Let  feL  [0,27 r] .  The  has  a  solution. 

* 

In  general  the  solution  of  Pf  is  not  uniaue.  Let  $  (f) 

*  * 

be  the  set  of  all  solutions  ( <j>  ,A  )  of  Pf.  Then  the  compactness 
of  U  and  the  continuity  of  T  implies  the  following  result  on 
continuous  dependence: 

2  *  *  * 

Theorem  6:  If  f_-*>f  in  L  [0,2tt],  (<£_, A_)e$  (f_),  then  there 
1  -  —  ■  n  n  n  n 

*  ★ 

exists  a  convergent  subsequence  of  and  everv  limit  point 

lies  in  (f ) . 

* 

Note  that  if  P^  is  uniquely  solvable,  i.e.  $  (f)  is  a 
*  *  , 

single  ordered  pair  C 4>  ,X  )  ,  then  Theorem  6  simply  says  that 

It  it  2 

f -*■  ( 4>  ,A  )  is  a  continuous  mapping  from  L  [0,2tt]  into  C(3D)xU. 

In  particular  the  operator  T  :f->-(<}>  ,A  ),  fcTU^X  satisfies  the 
~_1 

conditions  on  T  as  set  forth  in  the  opening  paragraph  of  this 

section.  If  Pf  is  not  uniquely  solvable  then  we  must  interpret 
T  as  a  set  valued  mapping  and  the  criteria  on  T  as  set  out 
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in  this  opening  paragraph  must  be  modified  in  an  obvious  manner. 

For  the  extension  of  the  above  results  to  the  case  of  the 
inverse  scattering  problem  for  electromagnetic  waves  see  [12]. 

The  inverse  scattering  problem  of  determining  the  shape 
of  a  "soft"  scattering  obstacle  from  a  knowledge  of  the  far  field 
pattern  can  be  handled  in  a  similar  manner  as  the  inverse  impedance 
problem  provided  we  can  define  an  appropriate  compact  family  of 
surfaces  U.  This  problem  was  discussed  in  [7]  where  the  set  U 
consisted  of  those  domains  D  contained  in  a  fixed  circular  annulus 
such  that  3D  has  a  uniformly  bounded  Holder  continuous  tangent. 

In  the  case  of  the  three  dimensional  inverse  scattering  problem 
it  was  also  necessary  to  assume  that  D  was  starlike  with  respect 
to  the  origin  ([1]). 

In  closing  we  mention  that  although  we  have  succeeded  in 
using  a  priori  information  to  stabilize  the  inverse  scattering 
problem,  we  have  made  no  statement  on  the  type  of  continuity 
that  results  (i.e.  Lipschitz,  Holder,  logarithmic,  etc.).  This 
is  an  important  problem  that  remains  to  be  investigated  since  if 
the  continuity  is  too  weak  there  is  little  hope  of  using  the 
optimization  problem  in  a  constructive  fashion  to  actually  deter¬ 
mine  X  or  D.  For  a  discussion  of  this  point  in  the  case  of 
linear  inverse  scattering  problems  we  refer  the  reader  to  [2]. 

—1 

V.  Constructive  Methods  for  Determining  l£x,xeTU. 

The  previous  section  shows  that  the  inverse  scattering 
problems  we  are  considering  in  this  paper  can  be  reformulated  as 
constrained,  nonlinear  optimization  problems.  The  numerical 
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solution  of  these  optimization  problems  is  presently  being 
investigated  by  A.  Kirsch  ([16]).  Related  approaches  have  been 
previously  considered  by  A.  Roger  ([18])  and  B.  Sleeman  ([19]). 
Needless  to  say  there  are  many  open  problems  connected  with  such 
a  numerical  approach,  particularly  in  connection  with  the  questions 
of  stability  and  convergence. 

The  computational  problems  associated  with  computing 

~_1 

Tx  ,xeTU,  are  considerably  simplified  if  one  has  access  to  low 
frequency  data.  We  first  briefly  consider  the  problem  of  deter¬ 
mining  the  surface  impedance  of  an  obstacle  from  low  frequency 
far  field  data  ([6],  [21]).  Let  F  be  the  far  field  pattern 
arising  from  the  addition  cf  the  scattered  waves  corresponding 
to  two  incoming  plane  waves  striking  the  obstacle  from  opposite 
directions.  If  we  expand  F  in  its  Fourier  series 

F ( 0 ;k)  *  ?  a  (k)ein9  (5.1) 

n=-°° 

then  the  (weighted)  low  frequency  limit  of  the  coefficient  aR(k) 

determine  a  secruence  of  numbers  b  from  which  we  can  define  the 

n 

harmonic  function 


h(x) 


00 

n=-°o 


V 


-  n 


.ine 


(5.2) 


It  can  be  shown  that  the  series  (5.2)  converges  for  r>_a  where 
a  is  the  radius  of  a  disk  containing  the  scattering  obstacle 
in  its  interior.  From  [6]  we  now  have  that  the  unknown  impedance 
X  is  the  solution  of  the  integral  ecruation  of  the  first  kind 


h(x) 


j  A  (y)  N  (x  ,y)  ds  (y) 
3D 


(5.3) 


.6 


where  N  denotes  the  Neumann  function  for  Laplace's  equation  in 
the  exterior  of  D.  Since  F  is  assumed  to  be  only  approximately 
known,  the  computation  of  the  coefficients  bn  (|n|<N  for  some 
positive  integer  N)  must  be  done  through  the  use  of  regularization 
procedures.  In  particular  the  left  hand  side  of  (5.3)  is  not 
known  exactly  and  the  problem  of  solving  (5.3)  is  hence  improperly 
posed.  However  stabilization  can  be  achieved  by  assuming  a  priori 
that  XeU ,  in  agreement  with  the  results  of  the  previous  section. 

A  similar  precedure  can  be  used  to  determine  the  shape  of 
a  soft  scattering  obstacle  from  a  knowledge  of  the  far  field 
pattern  at  low  values  of  the  frequency  ([3],  [9]).  This  method 
is  based  on  relating  the  low  frequency  limit  of  the  Fourier  coeffi' 
cients  of  the  far  field  pattern  to  the  coefficients  of  the  Laurent 
expansion  of  the  (unique)  analytic  function  f  that  conformally 
maps  the  exterior  of  the  unit  disk  onto  the  (unknown)  scattering 
obstacle  such  that  at  infinity  f  has  the  Laurent  expansion 

a.  a.~ 

f(w)  ■  aw  +  a  +  — —  +  +...  (5.4) 

o  w  z 

w 

where  a>0  is  the  mapping  radius.  In  particular  if  yn  denotes 
the  (weighted)  low  frequency  limit  of  the  Fourier  coefficients 
then  the  yn  are  related  to  the  an  by  a  relation  of  the  form 

u  »  2trn  an  ^  a  .  +  lower  order  coefficients,  (5.5) 

n  n-l 

for  n*l,2.,,,.  where  the  mapping  radius  a  can  be  determined  from 

experimental  data.  Furthermore,  if  the  coefficients  aQ,a^,...aN 

are  determined  from  (5.5)  and  f„  is  defined  bv 

N 


17 


f„(w)  =  aw  +  a  +  —  +  . 
N  o  w 


(5.6) 


then  from  the  Area  Theorem  in  geometric  function  theory  it  is 

2 

possible  to  deduce  the  L  -error  estimate 


<•*♦>  |2d*  <  jU-. 

o 


(5.7) 


Note  that  from  (5.5)  small  values  of  a  can  cause  large  errors  in 
the  computation  of  the  coefficients  a^  whereas  large  values  of 
a  imply  that  the  error  estimate  (5.7)  is  large.  Hence  in  order 
to  achieve  stability  a  must  be  bounded  from  above  and  below.  This 
will  be  guaranteed  if  D  is  known  a  priori  to  be  contained  in  a 
given  annulus,  in  agreement  with  the  results  of  the  previous  sec¬ 
tion.  Note  however  that  the  error  estimate  in  this  case  is  with 
2 

respect  to  the  L  norm  instead  of  the  pointwise  estimates  of 
Section  IV  and  hence  it  is  not  necessary  to  assume  a  priori  that 
3D  has  a  uniformly  bounded  Hdlder  continuous  tangent.  This  suggests 
the  problem  of  examining  the  stability  of  the  inverse  scattering 
problem  where  continuity  is  measured  with  respect  to  norms  different 
from  the  maximum  norm  and  determining  the  appropriate  compact  sets 
that  are  associated  with  these  norms. 


The  above  function-theoretic  approach  to  the  inverse 
scattering  problem  of  determining  the  shape  of  the  scattering 
obstacle  from  a  knowledge  of  the  far  field  pattern  was  first  given 
by  Colton  ([4])  and  Colton  and  Kleinman  ([9]).  Further  develop¬ 
ments  of  this  method  have  subsequently  been  provided  by  Hariharan 
([14]),  Sleeman  ([20]),  and  Smith  ([21]).  A  partially  successful 
attempt  to  extend  this  approach  to  the  inverse  scattering  problem 
in  l3  has  been  given  by  Colton  and  Kress  ([11]).  The  problem  in 
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this  case  is  that  instead  of  conformal  mapping  methods  one  must 
rely  on  the  use  of  level  curves  in  potential  theory  and  such  an 
approach  yields  much  weaker  results. 

The  discussion  provided  in  this  section  obviously  only 
touches  the  surface  of  the  computational  problems  inherent  in 
constructing  approximate  solutions  to  the  inverse  scattering 
problem  and  it  is  hoped  that  the  future  will  yield  new  developments 
and  insights  in  this  direction.  Indeed  the  effective  numerical 
solution  of  the  inverse  scattering  problem  is  the  basic  open 
problem  in  the  field  and  is  intimately  linked  to  all  of  the  pro¬ 
blems  we  have  previously  mentioned.  Our  hope  is  that  this  paper 
has  indicated  recent  progress  and  possible  new  directions  to  a 
sufficient  extent  to  encourage  others  to  enter  this  important  and 
fascinating  area  of  research. 
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